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Abstract 

We investigate the resummation effects in the Standard Model Higgs boson pair production at 
the LHC with soft-collinear effective theory. We calculate the invariant mass distribution and 
the total cross section at Next-to-Next-to-Leading-Logarithmic level with 7r 2 -enhanced terms re- 
summed, which are matched to the QCD Next-to-Leading Order results. Our results show that 
the resummation effects increase the Next-to-Leading Order results by about 30% ~ 40%, and the 
factorization scale uncertainty is reduced to 2% ~ 3%, which lead to increased confidence on the 
theoretical predictions. We also study the sensitivity of the total cross section and the invariant 
mass distribution to the Higgs boson self-coupling. We find that the total cross section and the in- 
variant mass distribution shape depend strongly on the Higgs boson self-coupling, and it is possible 
to extract Higgs boson self-coupling from the total cross section and invariant mass distribution 
when the measurement precision increases at the LHC. 
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I. INTRODUCTION 

Recently, at the CERN Large Hadron Collider (LHC), both the ATLAS [lj and CMS 
collaborations have found an about 5cr excess for a Standard Model (SM) Higgs boson 
like particle with a mass around 125 GeV. In the future, it is promising that LHC can 
tell whether this particle is the SM Higgs boson or something else when the measurement 
precision increases and its properties are probed. 

In the SM, the Higgs boson is responsible for the origin of electroweak symmetry breaking 
and the generation of elementary particle masses. After the Higgs field $ gets the vacuum 
expectation value v, the SM Higgs potential in the unitary gauge can be written as 



V(h) = A 



+ h) 2 v 



2 



(1) 



2 2 

where the Higgs boson self-coupling A is given by Asm — m 2 H /(2v 2 ) at the tree- level in the 
SM, and the radiative corrections can decrease Asm by 10% for ms = 125 GeV where main 
contributions from top quark loops [3] . In some new physics models, if a new heavy particle 
has the similar non-decoupling property to the top quark, much larger radiative corrections 
may exit. For example, in Two-Higgs-Doublet-Model, due to the non-decoupling effects of 
the additional heavier Higgs boson in loops, one-loop corrections can cause the lightest Higgs 
boson self-coupling to deviate from SM prediction by about 100%, even if the lightest Higgs 
boson couplings with gauge bosons and fermions, respectively, are almost SM-like [3] . Thus, 
it is important to measure the value of Higgs boson self-coupling to distinguish the SM from 
other models. 

At the LHC, the Higgs boson self-coupling A can be directly probed through Higgs boson 



pair production via gluon gluon fusion, and the relevant studies have been performed |4j-|l4|. 
In Ref. jsl, the Next-to-Leading Order (NLO) QCD corrections have been calculated in the 
large top quark mass limit, and the K-factor is found to be large (about 2), but the scale 
uncertainty at the NLO is still very high (about 30%). Hence, it is important to investigate 
higher order effects and perform QCD resummation calculations in order to improve the 
predictions on the cross section and reduce the theoretical uncertainties. Moreover, as is 
shown in Ref. 12] , the Higgs boson pair invariant mass distribution strongly depends on the 
Higgs boson self-coupling A. Thus in order to improve the precision of theoretical predictions 
it is also necessary to calculate QCD resummation effects on the invariant mass distribution. 
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Besides, compared to single Higgs boson production, due to the larger invariant mass in 
the final state, gluons radiated from initial states are much softer. Thus, the soft gluon 
resummation effects are expected to be much more important than in the case of single 
Higgs boson production. 

The full NLO QCD calculations for Higgs boson pair production including finite top 
quark mass effects are difficult. The approximated method of NLO QCD calculations ap- 



plied in previous works 



- |l4J | is multiplying the Leading-Order (LO) cross section for a 



finite top quark mass by a K-factor obtained in the large top quark mass limit. This ap 



15| , and we also perform 



proximated method work well in the single Higgs boson production 
resummation calculations in the large top quark mass limit. 

In this work, we present the threshold resummation in Higgs boson pair production at the 



LHC using soft-collinear effective theory (SCET) 16H20I]. which is very suitable to deal with 
the scattering processes with multiple scales. In the time-like process the threshold region is 
usually defined as the limit z = M 2 /s — > 1, where M is the invariant mass of the time-like 



particle and s is the square of the partonic center-of-mass energy 2lH29| . Generally, in the 
threshold region, the cross section can be factorized as 

a = U®S® f Pa ® f Pb , (2) 

where H, S and fp are the hard function, soft function and parton distribution func- 
tion (PDF), respectively. The hard function incorporates the short distance contributions 
arising from virtual corrections. The soft gluon effects coming from all colored particles 
are contained in the soft function. The PDF denotes the probability of finding a particular 
parton in the proton. In the SCET, the soft and collinear degrees of freedom can decouple 
via field redefinition, so the different scale physics can be separated [19| . One can deal with 
the relevant degrees of freedom at each scale and their scale dependences are controlled by 
remonetization group (RG) equations, respectively. As a result, the large logarithmic terms 
in the hard and soft functions can be resummed conveniently. 

The arrangement of this paper is as follows. In SecJIIl we give the fixed order results of 
Higgs boson pair production at the LHC. In Sec JIII|. we calculate the hard and soft matching 
coefficients at the NLO, and present RG-improved differential cross section analytically. In 
Sec JIVl we discuss the choice of hard and soft matching scale, and the numerical results of 
total cross section and the invariant mass distribution of Higgs boson pair. We conclude in 
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SecJVl 



II. FIXED ORDER RESULTS 




FIG. 1: Feynman diagrams for Higgs boson pair production at the tree level. 



In our calculations, we employ effective operators describing ggh and gghh interactions {f| 
where top quark has been integrated out. The relevant effective Lagrangian can be written 

as 

= ^C t (^)G^G« *>h - ^C t ^)G^G a ^h 2 (3) 
where the Higgs vacuum expectation value v is related to the Fermi constant by v — 

— 1/2 

{y2Gp\ . At the LO, Feynman diagrams for Higgs boson pair production induced by 
above effective Lagrangians are shown in Fig. [TJ Up to NLO, the Wilson coefficient C t (/J, 2 ) 
is not independent of the top quark mass, and can be written as 

«o,») = i + ^Q. (4) 

4 7T 

We consider Higgs boson pair production at the LHC with the center-of-mass energy Vs. 
The double differential cross section in the variables M 2 and Y, where M and Y respectively 



denotes the invariant mass and the rapidity of Higgs boson pair, can be written as 



30] 



d 2 a a 2 G 2 F M 2 I 4m| f 1 f 1 dz . ( ( , _ ( t/f . fK . 

dMW = 1152(2 7 r)3g V 1 ~ J * V J m T fi ^^ B ^^ Ci ^ V > M) (5) 

where fi/A(B) is the PDF for initial hadron A(B) and /// is the factorization scale. The hard 
subprocess kernel Cij(z,y, fif, M) depends on variable z and y, which are defined as 



M 2 Zl e -2Y _ z 



s (i_^)(i + a e - 



-2Y\ 



(6) 
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At the LO (ij) = (gg), and at the NLO (ij) = (gg), (qg), (qg), (gq), (gq), (qq), (qq). The 
parton momentum fractions x\ and x 2 can be written as 



r l-(l-y)(l-z) y It \-y(\-z) _ Y 
xi = \ ; — e , x 2 = \ -z ~ — re 



z l-y(l-z) 



zl-(l-y)(l-z) 



(7) 



In order to make sure that they are not exceed 1, the lower limit of integration of the variable 
z is defined as 



z(y) = max 



(re 2y - 1)(1 -y) + ^Ary 2 e 2Y + (1 - yf{\ - re 2Y ) 2 
(re- 2Y - l)y + yjgj - yfe~ 2Y + y 2 (l - re~ 2y ) 2 ' 

2(1 -y) 



(8) 



where the variable r is defined as r = M 2 / S. 




FIG. 2: Feynman diagrams for Higgs boson pair production at the NLO. 



Up to the NLO, the explicit results of Cij(z, y,fJ>f, M) can be written as follows, for which 



the corresponding Feynman diagrams are shown in Fig. [2j 
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(10) 
(11) 



M 2 — m 2 H + iYnmH 

where we set the renormalization scale to be the factorization scale fif. tuh and are the 
mass and total decay width of the Higgs boson, respectively. Here D(z) denotes the plus 
distribution as follows, 



D(z) 



1 , M 2 (l-2) 2 
in - 



1-z 



fJ? f Z 



(12) 



The other Cm can be obtained from the symmetry relations 



an ^aa-, ^ an ^ an a 



Con — Cn 



(13) 



We can split the integral kernel C gg as a sum of singular terms and regularized terms [22] 



C gg {z,y,iif, M) = + y) C(z,^,M) 
+C r ^(z,y, H ,M), 



6\v 2 



M 2 — m 2 H + iThtuh 



(14) 



where the singular terms consist of the contributions from the virtual and real soft gluons, 
and have the explicit form 

11 



C(z,(i f ,M) = 5(l-z) + — 

IT 



2tt 2 + — 5{1 - z) + 6D(z) 



(15) 
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After performing the integration over y, we can get the differential cross section in the 
threshold region 



2 ^.thresh 



d 2 a 



a 2 (3G 2 F M 2 

dM 2 dY ~ 1152(2vr) 3 5 



6Xv 2 



M 2 — m 2 H + iYhttih 



[ ^f 9 /A(^fe Y ,fi f )f g/B (V^e- 



Y 



/ —f a /A{Vre Y ,fif)f g/B (—e Y ,n f ) 

(16) 



y/re 

C(z,n f ,M) 



It is important to compare the contributions from the singular terms and NLO results. In 
Fig. [3]we see that the singular terms contribute 98% of the NLO total cross section and 95% 
at M = 500 GeV in Higgs boson pair invariant mass distribution. Thus the NLO results 
are well approximated by the singular terms, which reveals the necessity of resuming these 
singular terms to all orders. In Fig. |3]we also present our value of the total cross section (red 
solid line) and the results given by Eq. (21) in Ref. p] (black dashed line) at QCD NLO level. 
They are well consistent in the range of Monte Carlo integration error. 



_l 1 1 1 1 1 1 1 1 1 1 1 1 1 






MSTW2008nnlo ! 




iS=U TeV E 




-lo : 




NLO - 




singular terms _z 




S.Dawson et. al. J 


: , , 1 , , , , 1 


, , , , i , , , , i , , , , i , , , , i , , , , i , , , , i , , , ,- 



MSTW2008nnlo 
\fS=14 TeV 
- - LO 

NLO 

singular terms 



20 121 122 123 124 125 126 127 128 129 
m H (GeV) 




400 600 800 1000 1200 1400 1600 1800 2000 
M (GeV) 



FIG. 3: Comparison of the NLO results and singular terms contributions to the total cross section 
(Left) and the invariant mass distributions (Right) for Higgs boson pair production at the LHC 
with ^f~S = 14 TeV. The PDF is MSTW2008nnlo, the factorization scale is chosen as Higgs boson 
pair invariant mass. The other SM input parameters are shown in Eq. $4 



III. THE HARD AND SOFT FUNCTION AT NLO 



Following the approach in Ref. [24j, the hard scattering kernel C(z, fif, M) can be fac- 
torized as 



C{z,ti f ,M) = C 2 (n})n(M 2 ^})S - z) 2 ^ 2 ) 



(17) 



where Ct is the Wilson coefficient in Eq.(j4]), and it satisfies the RG equation 24 j 

d C t (^) = 1 t (a s )C t ( f i 2 ), with 7*(a.) = « s 2 ^-%^, (18) 



am/x aa s a 

where P(a s ) is the QCD /3-function. Then we obtain the evolution equation as 



^feSfeS^ (19) 



where /x t is the matching scale. "H(M 2 ,/x 2 ) and <S (s(l — z) 2 ,/x 2 ) in Eq. ( TT7|) are the hard 
and soft function, respectively. The hard and soft functions describe interactions at different 
scales, and they can be calculated order by order in perturbative theory at each scale. At 
the next-to-next-to-leading logarithmic (NNLL) accuracy, we need the explicit expressions 
of hard and soft functions up to NLO. 



Hard function 



In this process the hard function 7i(M , (jl . ) is the absolute value squared of the Wilson 
coefficient Cs{—M 2 ,fi 2 ) (here and below the negative arguments are understood with a — ie 



prescription) 



U{M*,n i ) = \C s 



-M 2 ,/i 2 )| 2 



(20) 



and Cs{—M 2 ,fi 2 ) can be obtained by matching the two-gluon operators in the full theory 
onto an operator in SCET, where the infrared divergences are subtracted in the MS scheme. 
We get the Wilson coefficient C 5 (— M 2 ,/i 2 ) up to NLO as follows, 



C s (-M 2 ^ 2 ) = l + 



« s (/i 2 
Air 



7T 



(21) 



where L = ln(— M 2 /fi 2 ), and this result agrees with those in Ref. 24|. The RG equation 
for Cs{—M 2 ,fi 2 ) is governed by the anomalous-dimension, the structure of which has been 
predicted up to four-loop level for the case of massless partons [31] . The Cs(—M 2 ,fi 2 ) 
satisfies the RG equation 



d 



dhxji 



C s (-M 2 ^ 2 ) 



cusp 



[a s ) In 



-M 2 
1^ 



+ 7 5 («s) 



C s (-M 2 ^ 2 ) 



(22) 



where 
ments 



(a s ) is the cusp anomalous dimension of Wilson loops with light-like seg- 
32], while 7 S (a s ) control the single-logarithmic evolution. After solving the RG 



equations, we get the Wilson coefficient Cs 



C s (-M 2 ,fi 2 f )=exp 



-M 2 

2S(/4,/iJ) - a r (/4,/iJ)ln — 2 S s (/4/4) 



C s (-M 2 ,fi 2 h ), (23) 



where fit is the hard matching scale, and 5'(z/ 2 ,/i 2 ) and ar(^ 2 ,/i 2 ) are defined as 



^ 2 ,/i 2 



a r (z/ 2 ,/x 2 ) 



a s (p 2 ) rA 



[a] 

a /?(«) 



(24) 
(25) 



a 7 s has a similar expression. Up to NNLL, we need 3-loop cusp anomalous dimension 
and 2-loop normal anomalous dimension, and their explicit expressions are collected in the 
Appendix of Refs. 



22 



24]. 



In Eq. fTST]) . if we set the hard matching scale /x 2 ~ M 2 , the perturbative behavior of fixer- 
order hard matching coefficient is bad due to the existence of 7r 2 terms from the squared 
logarithmic term. This poor perturbative behavior can be avoided if we set the hard match- 
ing scale at the space-like region fi\ ~ — M 2 . In order to get the hard matching coefficient 
at the time-like region, we apply the solution of RG Eq. ( 1221) to evolve the hard match- 
ing coefficient from the space- like region to time-like region. As a result, the 7r 2 — enhanced 
terms are resummed to all orders [23j. In the process of calculation we need to apply the 
strong coupling a s (u 2 ) evaluated at the negative arguments, where the explicit repressions 
are shown in Ref. 231 ] . 



B. Soft function 



The soft function S(s(l—z) 2 , /i 2 ), describing soft interactions between all colored particles, 



can be defined as[24j 



S(s(l - z) 2 , /i 2 ) = V~sW(s(l - z) 2 , /i 2 ). (26) 

In order to derive the RG equation satisfied by the soft function, we employ the fact that 
the physical cross section in the threshold region is independent of the arbitrary scale \i. 
We arrive at the integro-differential evolution equation for the momentum-space Wilson 
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coefficient 



d 



W(uj 2 ,^ 



din// 



[4T cnsp (a s )\n^- + 2 1 w (a s )}W(u 2 ,^ 

dJ W^' 2 ^ 2 )-W(co 2 ^ 2 ) 



^4r cuS p (<y s 



CO — U)' 



where the anomalous dimensions is 



(27) 



(28) 



in which 7 s is the coefficient of 6(1 — z) in the Altaralli-Parisi splitting function P gg (z) 22]. 
The solution of this evolution equation is 24] 



CO 



W(u 2 ,fi 2 f ) = exp [-4£(nl,n 2 f ) + 2a 7 w(/i 2 , $)] s(d v ,/i 2 ) ( — 



2 \ V p -2j EV 



T(2 V ) 



(29) 



where 77 = 2ar(//g,/^), and s(L s ,fi 2 ) is the Laplace transformation of the soft Wilson loop 
at the matching scale y, s . At the NLO it is given by 



KL„tf) = 1 + ^11(6^ + tt 2 ). 



(30) 



Eq. ( 129]) is well defined when 77 > 0. When 77 < the results should be understood with the 
pole contribution subtracted, and we deal with it following the approach in Ref. 22]. 



C. Final RG improved differential cross section 



After combining the hard and soft function together, we obtain the RG-improved expres- 
sion for Eq.( lT5|) . 



xs In 



2 -?)] 2 
M 2 (l-z) 2 



/j? s z 



(31) 



where 



a 2.2 2 2, a 2 (fi 2 ) \(3(a s (n 2 ))/a 2 (fi 2 ' 



— M 5 



2 \ -2ar(l*£,J*«) 



a 2 s (n 2 ) [(3(a s (fi 2 t ))/a 2 s (fi 2 ) 
|exp [AS(n 2 h ,fi 2 ) - 2a yS (fi 2 h , fi 2 ) + Aa yB (fi 2 s , fi 2 f )] \ . (32) 



In the above expression, if we set scales fi t = [i>h = A*s — A 4 /, then we can recover the 
threshold singular terms in Eq. ffT5]) . which appear in the fixed-order calculation. Substituting 
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Eq.( l3T|) into Eq.( IT6|) . we can get the resummed differential cross section at the threshold 
region. 

To give precise predictions, we resum the singular terms to all orders and include the 
non-singular terms up to NLO. The RG-improved differential cross section can be written 

as 



d 2 a Resum _ d 2 a thTCS ( d 2 (x NLO dV hres \ 
dM 2 dY ~ dM 2 dY + [ dM 2 dY ~ dM 2 dY ) 



(33) 

expanded to NLO 



In the threshold region, the second term in the above expression vanishes, and only the 
threshold contribution dominates, and in the region far away from the threshold region, the 
resummation effect is not important and the differential cross section is dominated by the 
NLO results. 



IV. NUMERIAL DISCUSSION 



In this section, we discuss the numerical results for the threshold resummation effects on 
Higgs boson pair production at the LHC We choose the following SM imput parameters 

G F = 1.16637 x 1(T 5 GeV -2 , m H = 125 GeV, T H = 4.03 x 1(T 3 GeV, A = A SM - (34) 

Throughout the numerical calculations, we use the MSTW2008nnlo PDF sets and associated 
strong coupling constant. The factorization scale /i/ is choose as the Higgs boson pair 
invariant mass M unless specified otherwise. There are still three matching scales, i.e. fit, 
fih, ^si which should be choose such that the Wilson coefficients evaluated at the matching 
scales have stable perturbative expansions. 



A. Scale setting 



Because the effective operator is constructed when top quark loop is integrated out, it is 



natural to choose the matching scale as fi t = Tn t 24j 
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-\i 2 /M 2 



FIG. 4: Dependence of c\ on the hard matching scale squared 



The hard matching coefficient can be expanded in the perturbative series as 

C s {-M\ /4) = 1 + ^ cj ^ <«■ (35) 

n=l ^ ' 

In Fig. H] the dependence of c\ on the hard matching scale squared is shown. The 
coefficient c\ vanishes for —y^\ ~ 3.6M 2 and — ~ 0.28Af 2 . Around the region of second 
solution the coefficient C\ varies strongly, thus we discard it. In our numerical results, we will 
take — ix\ = 3.6M 2 as the default choice. Thus the ir 2 — enhanced corrections are contained 
into the evolution function U. 

For the soft matching scale fi s , the situations are not so clear, since the soft function 
<S(s(l — z) 2 ,/i 2 ) depends on the variable z. The soft scale is chosen such that the pertur- 
bative expression of the soft function have a well-behaved convergence after performing the 
integration over the variable z. 
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FIG. 5: Left: The contributions from soft functions to the NLO cross section. Right: Results for 
the soft matching scale fj, s for different r. 

In Fig [51 we present the corrections from soft functions at the NLO. Following the methods 
in Refs. |22|, [24|, we determine the soft matching scale as follows. 

• Starting from a high scale, choose the value of // s where the soft correction drops below 
40%; 

• Choose the value of fi s where the soft correction is minimized. 
As a result, we get two soft scales: 

l + 60r ' " V7 + 400F { } 

In Fig. El we show the dependence of the resummation cross section at the LHC with 
y/S = 14TeV on the hard matching scale and the soft matching scale, where the two match- 
ing scales are changed in the region M 2 < —fi 2 < 9M 2 and fi 1 / < fi s < /if, respectively. 
From Fig. [61 it is obvious that the scale dependence is reduced at the NNLL, and we can 
neglect their impacts on the cross section. 
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FIG. 6: The scale dependence of the Higgs boson pair production resummation cross section at 
the NLL and NNLL on the hard matching scale (Left) and the soft matching scale (Right). 

The factorization scale dependence of the total cross section is shown in Fig. [7J where 
the factorization scale fif is changed from M/2 to 2M. In the fixed-order calculation, 
the scale dependence is very large, and even at the NLO the scale uncertainty is about 
30%. However, the factorization scale dependence is obviously reduced after performing 
resummation calculations, and the scale uncertainty is only 2 ~ 3% in the NNLL. 




FIG. 7: The factorization scale dependence of the Higgs boson pair production fixed-order (Left) 
and resummed (Right) total cross section . 
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B. Higgs boson pair invariant mass distributions 



0J5 i i i i i i i i i i i i i i i i i i i i i i i i i i i ' i ' 




M (GeV) M (GeV) 



FIG. 8: Invariant-mass distributions for Higgs boson pair production at the LHC with \/S = 14 TeV 
and 33 TeV. The bands indicate the factorization scale uncertainty. The blue and red bands are 
fixe-order results at LO and NLO, respectively, and the green band includes the effects of NNLL 
resummation matched to NLO results. 

Although the calculations in the large top quark mass limit are not sufficient to reproduce 
all kinematics distribution of the full theory, it is also significant to study the enhancement 
effects of resummation on the kinematics distribution. In Fig. [HJ we present the Higgs boson 
pair invariant mass distribution at the LHC with \/S = 14 and 33 TeV, where the band 
reflects the theoretical uncertainty by varying the factorization scales between M/2 and 
2M. It is obvious that the resummation calculations increase the differential cross section 
da/dM 2 by 30% at M = 500 GeV and 40% at M = 1500 GeV, respectively, with VS = 14 
TeV. Moreover, the factorization scale uncertainties are reduced to about 3% at M — 500 
GeV and in the large invariant mass region the scale uncertainties almost vanish. 

In order to estimate the effects of soft-gluon resummation, we define the K-factor of the 
invariant mass distributions as 



da 



dM 2 



K(M)-^- ■ (37) 

NLO 



Resum dM 2 

In Fig. we present the K-factor for different Higgs boson pair invariant mass, where the 
bands also reflect the scale uncertainties. Obviously, the K-factors increase with the increas- 
ing of the invariant mass, because gluons radiated from initial states are softer in the large 
invariant mass region than in the small invariant mass region. 
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FIG. 9: K-factor for Higgs boson pair invariant mass distribution at the LHC with y/S = 14 TeV 
and 33 TeV, and the bands indicate the factorization scale uncertainty. 



C. Total cross section 

After performing integration over the Higgs boson pair invariant mass M, we can get the 
total cross sections. Table [T]shows the total cross section at the LHC with v^S* = 8, 14 and 33 
TeV. Our results show that the resummation effects increase the NLO results by about 
30% ~ 40% and the factorization scale uncertainty is reduced from 30% to 3%, which leads 
to increased confidence on the theoretical predictions. 
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TABLE I: Total cross sections of Higgs boson pair production at the LHC in fb for different Higgs 
mass. The errors are the factorization scale uncertainties. 
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D. Higgs boson self-coupling measurements at the LHC 




FIG. 10: The dependence of total cross section on the Higgs boson self-coupling A at the LHC. 
The bands indicate the factorization scale uncertainty. 

After the discovery of the Higgs-like particle is confirmed, it is important to verify whether 
it is the SM Higgs or not. Thus it is necessary to measure the Higgs boson self-coupling A at 
the LHC through Higgs boson pair production. In the expression of effective Lagrangian ([3]) 
the first and the second operator induce the first and the second Feynman diagram in Fig. HJ 
respectively, where only the physical process represented by the first Feynman diagram is 
sensitive to the value of Higgs boson self-coupling. Due to the different sign of these two terms 
in Eq. ([3]) important interference effects exist. If A < 0, the contributions from these two 
terms can be added, and if A > 0, the contributions cancel each other. Thus the total cross 
section of Higgs boson pair production decreases as increasing the value of A. This relation 
is shown in Fig. [TUl where we present the dependence of the Higgs boson pair production 
cross section on the self-coupling A. Obviously, the theoretical predictions at the NLO level 
are not sufficient to precisely determine the value of A due to high scale dependence and 
underestimate of the total cross section. After including the NNLL resummation effects, it 
is promising to distinguish the SM Higgs boson self-couplings from non-SM ones. 

Moreover, this important interference effects also contribute to the invariant mass dis- 
tribution of the Higgs boson pair production. In Fig. [HJ we present the normalized Higgs 
boson pair invariant mass distribution at the LHC with = 14 TeV where the scale un- 
certainties almost vanish due to the cancelation of scale uncertainties in ratio. When A < 0, 
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the peak of the distribution locates at the small invariant mass region (yellow solid line). 
The position of peak moves towards large invariant mass region with the increasing of the 
value of A. Especially, when A > Asm there exits two different peak (blue dashed line). Here, 
we explain the reason why the shape of the invariant mass distribution strongly depends on 
the value of A. As is shown in Eq.f ll6j) the differential cross sections share a common factor 

i - 6Xy2 (m 

M 2 -m 2 H + im H T H ' K ' 

where the second terms represent a Higgs boson exchanged in the s-channel. When A < 0, 
the invariant mass distributions at the small value region are added up, thus an obvious 
peak exists in the small invariant mass region. When A > 0, in the small invariant mass 
region the contributions from two terms cancel each other, so the peak of the invariant mass 
distribution moves toward large value region with the increasing value of A. In particular, 
when M 2 = 6Xv 2 + m 2 H the common factor equals zero. Considering the threshold limit 
M > 2m h, there exists a zero point in the invariant mass distribution as A > Asm- For 
example, as shown in Fig. [TTJ there are two peaks in the invariant mass distribution when 
A = 2 Asm, where the differential distribution equals zero with M = 330.7 GeV. From above 
discussions, we see that it is possible to extract the parameter A from the total cross section 
and Higgs boson pair invariant mass distribution when the measurement precision increases 
at the LHC. 



i 1 1 1 i 1 1 1 i 1 1 

NNLL+NLO 

A^ SM = 

-— ^Km = 1 

^SM = 2 




200 400 600 800 1000 1200 1400 1600 1800 2000 
M (GeV) 



FIG. 11: The normalized Higgs boson pair invariant mass distribution at the LHC with = 14 
TeV. 
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V. CONCLUSION 



We have calculated the resummation effects in the SM Higgs boson pair production at 
the LHC with SCET. We present the invariant mass distribution and the total cross section 
at NNLL with 7r 2 -enhanced terms resummed, which are matched to the NLO results. Our 
results show that the resummation effects increase the NLO results by about 30% ~ 40%, 
and the factorization scale uncertainty is reduced to 2% ~ 3%, which leads to increased 
confidence on the theoretical predictions. We also study the sensitivity of the total cross 
section and the invariant mass distribution to the Higgs boson self-coupling. We find that the 
NNLL resummation effects can decrease the needed luminosity to detect the Higgs boson 
self-coupling about 30%, and the total cross section and the invariant mass distribution 
shape depend strongly on the self-coupling. As a result, it is possible to extract Higgs boson 
self-coupling from the total cross section and Higgs boson pair invariant mass distribution 
when the measurement precision increases at the LHC. 
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